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ABSTRACT 

The Green's function for rf-driven currents in a 
toroidal plasma is calculated in the long mean-free- 
path Umit. 

INTRODUCTION 

The "adjoint" technique provides an elegant 
method for computing the current driven in a 
plasma by externally injected radio-frequency 
waves. This technique involves using the self- 
adjoint property of the linearized collision oper- 
ator to express the current in terms of a Green's 
function, where the Green's function is propor- 
tional to the perturbed distribution in the presence 
of an electric field. This technique was first ap- 
plied by Hirshman^ for the case of neutral-beam 
driven currents in a homogeneous magnetic field. 
The effect of the trapped particles found in toroidal 
confinement systems was included by Taguchi.^ 
The application of this technique to rf-driven cur- 
rents was made by Antonsen and Chu,"^ Taguchi,^ 
and Antonsen and Hui.^ Approximations to the 
Green's function have been given by a number of 
authors. 

In this paper, we solve the adjoint problem nu- 
merically in toroidal geometry using the collision 
operator for a relativistic plasma. ^^"^^ The per- 
tinent approximations are: the device is assumed 
to be axisymmetric; the mean-free-path is assumed 
to be long compared to the device (the "banana 
limit"); drifts of the electrons away from the ini- 
tial flux surface are neglected; in addition the ex- 
pansion of the collision operator in Legendre har- 
monics is truncated after the Pi (cos 9) term. By 
posing the problem in terms of a Green's function, 
we are, of course, also assuming that the plasma is 
close enough to equilibrium for the collision oper- 
ator to be linearized, and that the wave-driven flux 
is known. 

BASIC EQUATIONS 



est order, a function of the collisionless constants 
of motion. We choose to express the distribution in 
terms of the "midplane" coordinates^"* (mq, 6*0), the 
magnitude and direction with respect to the mag- 
netic field of the momentum per unit mass (hence- 
forth called just the momentum) at the position 
where the magnetic field is minimum. Measuring 
position on the flux surface by the length I along 
the field line from this point, the momentum {u, 9) 
at an arbitrary position is 



"0, 



81-0^9 



bs'm 9o, 



where b = B{l)/B{0). Particles with sin^ 6*0 > 
sin^ 6'tr = l/6max ^re trapped; other particles are 
passing. Assuming that the rf is sufficiently weak, 
the distribution satisfies 
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where fmC{f / fm) is the collision operator lin- 
earized about a Maxwellian /,n, 

n J v\\ 

is the bounce-averaging operator, ti, = J dl/v\\ is 
the bounce time, A = ti,vq cos 6*0 / L, and L is total 
length of the field line (from one intersection with 
the midplane to the next). Sq is the rf-induced flux 
in momentum space expressed in midplane coordi- 
nates. This is related to the local rf-induced flux S 
via 
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We should also include a term in eq. (1) which 
reflects the slow heating of the background elec- 
trons.'' However, this term does not contribute to 
the current carried by /. 

The power dissipated by the wave between two 
neighboring flux surfaces is 



In the long mean-free-path limit, the electron 
distribution on a particular flux surface is, to low- 



(2) 



where dV is the elemental volume between the two 
surfaces and 



Po = m d uo ASo • vq. 



The current density at the midplane is 
"^011 =Q rf^uoi'ocos^o/- 



(3) 
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At an arbitrary point the current density is J|| = 
6 Jo II . The total toroidal current flowing between 
two neighboring flux surfaces is 



where 



Q 



QdV 
Jdl/b 



Oh 
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dl B, 



C 



2nR B 

is the safety factor, R is the major radius, and 
is the toroidal magnetic field. 

Rather than determine Jqh directly by solving 
eq. (1), we consider the adjoint problem. 



(C(X)) 



Vq cos ^0 ^ 
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where 6 = 1 for passing particles and for trapped 
particles. This is the equation for the perturbed 
electron distribution is the presence of a toroidal 
loop voltage TL/Q. The rf-driven current density 
is then given by^ 
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Joii = I d Uo ASo • 



(7) 



We will express the current drive efficiency by the 
ratio r] = Joh/Pq. Another useful measure of ef- 
ficiency is in terms of the macroscopic variables I 
and W, namely 
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W 
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BOUNCE-AVERAGED COLLISION 
OPERATOR 

The linearized colUsion operator is made up of 
three terms 

C{X) = (C^/^(/mX, /m) + C^/^(/m, /mX) 

+ C'^\fmXji))/U, (8) 



electron-ion term is computed in the Lorentz limit 
(with rrii oo). It can be combined with the first 
term and bounce averaged to give 
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where D^u and F„ are the coefficients of energy 
diffusion and drag due to electron-electron colli- 
sions, and De0 is the pitch-angle scattering coeffi- 
cient due to colUsions with both electrons and ions. 
These are given by one-dimensional integrals over 
a Maxwellian distribution.^'^ 

Let us now turn to the second term in eq. (8). 
Since x is odd in U||, we can expand x{u, 0, 1) in 
terms of spherical harmonics as foUows: 

x{u,e,i)= ^ xkiu,i)Pk{cose), 

k odd 

where Pk is the Legendre polynomial of degree k 
mdxk{u,l) = (2fc+l) j;/\{u,e,l)Pk{cose) X 
sin^rf^. The Unearized collision operator is a 
spherically symmetric, so that its angular eigen- 
functions are spherical harmonics. This allows us 
to vmte the term C^/'=(/m, Xk{u, Z)Pfe(cos6i))//m 
as Ck = Pk{cos 9)Ik{xk{u, 0)' where Ik is a lin- 
ear integral operator. Transforming to midplane co- 
ordinates, we find 



Xk{u,l) = {2k + 



rn/2 
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x{uo,6o)Pksm 6od0o, 



where Pk = Pfc(cos0)cos6'o/cos0 and we have 
used the fact that x is zero for trapped particles. 
The bounce-averaged colhsion term becomes 

{Ck) = \j ^PkhixkiuJ)). 

Evaluating these expressions is simplified by de- 
composing Pk into midplane Legendre harmonics: 

Pi = Pifi, 

P3 = -(6-l)Pi,o + 6P3,0, 



Pk 



Gk,k'Pk',0, 



fe' = l,3,... 



where C°'^''{fa, fb) is the collision operator for where P^ o = Pfe(cos^o) and Gk,k' is a polyno- 
distribution fa coUiding off distribution /;,. The mial in b. The colhsion term can now be written 
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as 
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where 



k',k" odd 



Hk,k'k" = r.7 ;/ , 1 l.'Gk,k" , 
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Xfe.o = Xk{l = 0), and j4 = J Adl/L. In particu- 
lar, we have 

(co = 6^/i(xi.o)e. 

At present, we include only the k — 1 term, ig- 
noring all terms {Ck>3)- We can estimate the error 
incurred by comparing the results we get for the 
electrical conductivity with those of Rosenbluth et 
al}^ This indicates that the relative error in x is on 
the order of O.OS-yi where e is the inverse aspect 
ratio. 

NUMERICAL SOLUTION 

Putting all the terms in eq. (6) together, we ob- 
tain for the passing particles 
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-cos^o r / X , t^ocos6io ^ 
+ ^^—^1 (Xi,o) + q — ^ — = 0. 
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We solve this integro-differential equation numeri- 
cally in the domain < Oq < 6ti, with boundary 
condition x{(^o = (^tr) = 0. 

A simple magnetic field configuration with cir- 
cular flux surfaces is chosen. Designating the 
poloidal angle by (p, we choose 

R = -Ro(l + ecos(p), 
Be =Sfo/(l + ecos<?i), 
= B^q/{1 + ecos(j)), 
b= (1 -|-e)/(l-|-ecos0). 

This gives l/L = (f)/2n, Q =_{B^o/Bj,q) x 

e/VT^, and L = 2n R^Q ^l + ^J^^. 

We normalize velocities and momenta to the 
thermal speed Ut = \JT jm, times to the in- 
verse collision frequency v'^ ^ , where vt = ng^ x 
log A/47reom^Uj, x to qutjvu efficiency j] to 



q/mutvt. The plasma is characterized three di- 
mensionless parameters: T/ mc?, e, and the effec- 
tive ion charge state Z. 

Level curves for x(^*Oj ^0) for ^ typical case are 
shown in fig. 1. In computing the efficiency we 
specialize to waves which push the particles paral- 
lel to the magnetic field. First, we suppose that the 
wave is absorbed in a single location in momentum 
space, i.e.. So cx (5(uo — Uq)uo||. The efficiency 
is given by = {dx/duQ\\)lvQ\\ evaluated at Uq. 
This is shown as a function of Ug in fig. 2. This 
shows where in momentum space we should try to 
have waves absorbed in order to maximize the effi- 
ciency. 

We consider current drive by waves which are 
Landau damped. We assume that the rays pierce 
the flux surface at a single poloidal angle (j)' where 
b = b', and that the wave does not alter the 
slope of the electron distribution appreciably. In 
this case, we have S oc fniS{v\\ — Uph)^(0 — 
^')u|| and the bounce-averaged flux is given by 
ASq oc fmS{vo COS0' — i;ph)uo||> where sin^ 9' = 
b' sin^ 60. The current drive efficiency can be cal- 
culated by inserting this form for So into eqs. (3) 
and (7) and performing the integrals numerically. 
The resulting efficiencies are given in fig. 3. This 
confirms that toroidal effects reduce the efficiency 
of current drive and that this reduction can be min- 
imized if the waves are absorbed on the high-field 
side of the torus where there are fewest trapped par- 
ticles. 
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Figure 2: Contour plot of efficiency rj of current 
drive with point excitation of the rf, i.e.. So oc 
5{uo - u'q)uo\\. Here, Z = 1, e = 0.03, and 
T/mc^ — 0.05. The levels of the contours are 
given by ry = 2j for integer j > 1 increasing from 
the origin outwards. 
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Figure 3: Efficiencies for current drive by Landau- 
damped waves for Z = 1, and Tjmc? = 0.05. 
The top curve gives the efficiency for the case of a 
uniform magnetic field e = 0.0. The other curves 
are for e = 0.03 and three different poloidal angles 
at which the wave is absorbed. 
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